Inferring COVID-19 testing and vaccination behavior from New Jersey testing data.

Characterizing the relationship between disease testing behaviors and infectious disease dynamics is of great importance for public health. Tests for both current and past infection can influence disease-related behaviors at the individual level, while population-level knowledge of an epidemic's course may feed back to affect one's likelihood of taking a test. The COVID-19 pandemic has generated testing data on an unprecedented scale for tests detecting both current infection (PCR, antigen) and past infection (serology); this opens the way to characterizing the complex relationship between testing behavior and infection dynamics. Leveraging a rich database of individualized COVID-19 testing histories in New Jersey, we analyze the behavioral relationships between PCR and serology tests, infection, and vaccination. We quantify interactions between individuals' test-taking tendencies and their past testing and infection histories, finding that PCR tests were disproportionately taken by people currently infected, and serology tests were disproportionately taken by people with past infection or vaccination. The effects of previous positive test results on testing behavior are less consistent, as individuals with past PCR positives were more likely to take subsequent PCR and serology tests at some periods of the epidemic time course and less likely at others. Lastly, we fit a model to the titer values collected from serology tests to infer vaccination trends, finding a marked decrease in vaccination rates among individuals who had previously received a positive PCR test. These results exemplify the utility of individualized testing histories in uncovering hidden behavioral variables affecting testing and vaccination.

Polarization and the Psychology of Collectives.

Achieving global sustainability in the face of climate change, pandemics, and other global systemic threats will require collective intelligence and collective action beyond what we are currently experiencing. Increasing polarization within nations and populist trends that undercut international cooperation make the problem even harder. Allegiance within groups is often strengthened because of conflict among groups, leading to a form of polarization termed "affective." Hope for addressing these global problems will require recognition of the commonality in threats facing all groups collective intelligence that integrates relevant inputs from all sources but fights misinformation and coordinated, cooperative collective action. Elinor Ostrom's notion of polycentric governance, involving centers of decision-making from the local to the global in a complex interacting framework, may provide a possible pathway to achieve these goals.

On the limits to invasion prediction using coexistence outcomes.

The dynamics of ecological communities in nature are typically characterized by probabilistic processes involving invasion dynamics. Because of technical challenges, however, the majority of theoretical and experimental studies have focused on coexistence dynamics. Therefore, it has become central to understand the extent to which coexistence outcomes can be used to predict analogous invasion outcomes relevant to systems in nature. Here, we study the limits to this predictability under a geometric and probabilistic Lotka-Volterra framework. We show that while individual survival probability in coexistence dynamics can be fairly closely translated into invader colonization probability in invasion dynamics, the translation is less precise between community persistence and community augmentation, and worse between exclusion probability and replacement probability. These results provide a guiding and testable theoretical framework regarding the translatability of outcomes between coexistence and invasion outcomes when communities are represented by Lotka-Volterra dynamics under environmental uncertainty.

Fundamental bound on epidemic overshoot in the SIR model.

We derive an exact upper bound on the epidemic overshoot for the Kermack-McKendrick SIR model. This maximal overshoot value of 0.2984 · · · occurs at [Formula: see text]. In considering the utility of the notion of overshoot, a rudimentary analysis of data from the first wave of the COVID-19 pandemic in Manaus, Brazil highlights the public health hazard posed by overshoot for epidemics with R0 near 2. Using the general analysis framework presented within, we then consider more complex SIR models that incorporate vaccination.

Coupled models of genomic surveillance and evolving pandemics with applications for timely public health interventions.

Disease surveillance systems provide early warnings of disease outbreaks before they become public health emergencies. However, pandemics containment would be challenging due to the complex immunity landscape created by multiple variants. Genomic surveillance is critical for detecting novel variants with diverse characteristics and importation/emergence times. Yet, a systematic study incorporating genomic monitoring, situation assessment, and intervention strategies is lacking in the literature. We formulate an integrated computational modeling framework to study a realistic course of action based on sequencing, analysis, and response. We study the effects of the second variant's importation time, its infectiousness advantage and, its cross-infection on the novel variant's detection time, and the resulting intervention scenarios to contain epidemics driven by two-variants dynamics. Our results illustrate the limitation in the intervention's effectiveness due to the variants' competing dynamics and provide the following insights: i) There is a set of importation times that yields the worst detection time for the second variant, which depends on the first variant's basic reproductive number; ii) When the second variant is imported relatively early with respect to the first variant, the cross-infection level does not impact the detection time of the second variant. We found that depending on the target metric, the best outcomes are attained under different interventions' regimes. Our results emphasize the importance of sustained enforcement of Non-Pharmaceutical Interventions on preventing epidemic resurgence due to importation/emergence of novel variants. We also discuss how our methods can be used to study when a novel variant emerges within a population.

Role of masks in mitigating viral spread on networks.

Masks have remained an important mitigation strategy in the fight against COVID-19 due to their ability to prevent the transmission of respiratory droplets between individuals. In this work, we provide a comprehensive quantitative analysis of the impact of mask-wearing. To this end, we propose a novel agent-based model of viral spread on networks where agents may either wear no mask or wear one of several types of masks with different properties (e.g., cloth or surgical). We derive analytical expressions for three key epidemiological quantities: The probability of emergence, the epidemic threshold, and the expected epidemic size. In particular, we show how the aforementioned quantities depend on the structure of the contact network, viral transmission dynamics, and the distribution of the different types of masks within the population. Through extensive simulations, we then investigate the impact of different allocations of masks within the population and tradeoffs between the outward efficiency and inward efficiency of the masks. Interestingly, we find that masks with high outward efficiency and low inward efficiency are most useful for controlling the spread in the early stages of an epidemic, while masks with high inward efficiency but low outward efficiency are most useful in reducing the size of an already large spread. Last, we study whether degree-based mask allocation is more effective in reducing the probability of epidemic as well as epidemic size compared to random allocation. The result echoes the previous findings that mitigation strategies should differ based on the stage of the spreading process, focusing on source control before the epidemic emerges and on self-protection after the emergence.

Medium-term scenarios of COVID-19 as a function of immune uncertainties and chronic disease.

As the SARS-CoV-2 trajectory continues, the longer-term immuno-epidemiology of COVID-19, the dynamics of Long COVID, and the impact of escape variants are important outstanding questions. We examine these remaining uncertainties with a simple modelling framework that accounts for multiple (antigenic) exposures via infection or vaccination. If immunity (to infection or Long COVID) accumulates rapidly with the valency of exposure, we find that infection levels and the burden of Long COVID are markedly reduced in the medium term. More pessimistic assumptions on host adaptive immune responses illustrate that the longer-term burden of COVID-19 may be elevated for years to come. However, we also find that these outcomes could be mitigated by the eventual introduction of a vaccine eliciting robust (i.e. durable, transmission-blocking and/or 'evolution-proof') immunity. Overall, our work stresses the wide range of future scenarios that still remain, the importance of collecting real-world epidemiological data to identify likely outcomes, and the crucial need for the development of a highly effective transmission-blocking, durable and broadly protective vaccine.

Epidemiological impacts of post-infection mortality.

Infectious diseases may cause some long-term damage to their host, leading to elevated mortality even after recovery. Mortality due to complications from so-called 'long COVID' is a stark illustration of this potential, but the impacts of such post-infection mortality (PIM) on epidemic dynamics are not known. Using an epidemiological model that incorporates PIM, we examine the importance of this effect. We find that in contrast to mortality during infection, PIM can induce epidemic cycling. The effect is due to interference between elevated mortality and reinfection through the previously infected susceptible pool. In particular, robust immunity (via decreased susceptibility to reinfection) reduces the likelihood of cycling; on the other hand, disease-induced mortality can interact with weak PIM to generate periodicity. In the absence of PIM, we prove that the unique endemic equilibrium is stable and therefore our key result is that PIM is an overlooked phenomenon that is likely to be destabilizing. Overall, given potentially widespread effects, our findings highlight the importance of characterizing heterogeneity in susceptibility (via both PIM and robustness of host immunity) for accurate epidemiological predictions. In particular, for diseases without robust immunity, such as SARS-CoV-2, PIM may underlie complex epidemiological dynamics especially in the context of seasonal forcing.

Individual costs and societal benefits of interventions during the COVID-19 pandemic.

Individual and societal reactions to an ongoing pandemic can lead to social dilemmas: In some cases, each individual is tempted to not follow an intervention, but for the whole society, it would be best if they did. Now that in most countries, the extent of regulations to reduce SARS-CoV-2 transmission is very small, interventions are driven by individual decision-making. Assuming that individuals act in their best own interest, we propose a framework in which this situation can be quantified, depending on the protection the intervention provides to a user and to others, the risk of getting infected, and the costs of the intervention. We discuss when a tension between individual and societal benefits arises and which parameter comparisons are important to distinguish between different regimes of intervention use.

Spreading processes with mutations over multilayer networks.

A key scientific challenge during the outbreak of novel infectious diseases is to predict how the course of the epidemic changes under countermeasures that limit interaction in the population. Most epidemiological models do not consider the role of mutations and heterogeneity in the type of contact events. However, pathogens have the capacity to mutate in response to changing environments, especially caused by the increase in population immunity to existing strains, and the emergence of new pathogen strains poses a continued threat to public health. Further, in the light of differing transmission risks in different congregate settings (e.g., schools and offices), different mitigation strategies may need to be adopted to control the spread of infection. We analyze a multilayer multistrain model by simultaneously accounting for i) pathways for mutations in the pathogen leading to the emergence of new pathogen strains, and ii) differing transmission risks in different settings, modeled as network layers. Assuming complete cross-immunity among strains, namely, recovery from any infection prevents infection with any other (an assumption that will need to be relaxed to deal with COVID-19 or influenza), we derive the key epidemiological parameters for the multilayer multistrain framework. We demonstrate that reductions to existing models that discount heterogeneity in either the strain or the network layers may lead to incorrect predictions. Our results highlight that the impact of imposing/lifting mitigation measures concerning different contact network layers (e.g., school closures or work-from-home policies) should be evaluated in connection with their effect on the likelihood of the emergence of new strains.

Challenges in cybersecurity: Lessons from biological defense systems.

Defending against novel, repeated, or unpredictable attacks, while avoiding attacks on the 'self', are the central problems of both mammalian immune systems and computer systems. Both systems have been studied in great detail, but with little exchange of information across the different disciplines. Here, we present a conceptual framework for structured comparisons across the fields of biological immunity and cybersecurity, by framing the context of defense, considering different (combinations of) defensive strategies, and evaluating defensive performance. Throughout this paper, we pose open questions for further exploration. We hope to spark the interdisciplinary discovery of general principles of optimal defense, which can be understood and applied in biological immunity, cybersecurity, and other defensive realms.

Evolutionary dynamics within and among competing groups.

Biological and social systems are structured at multiple scales, and the incentives of individuals who interact in a group may diverge from the collective incentive of the group as a whole. Mechanisms to resolve this tension are responsible for profound transitions in evolutionary history, including the origin of cellular life, multicellular life, and even societies. Here, we synthesize a growing literature that extends evolutionary game theory to describe multilevel evolutionary dynamics, using nested birth-death processes and partial differential equations to model natural selection acting on competition within and among groups of individuals. We analyze how mechanisms known to promote cooperation within a single group-including assortment, reciprocity, and population structure-alter evolutionary outcomes in the presence of competition among groups. We find that population structures most conducive to cooperation in multiscale systems can differ from those most conducive within a single group. Likewise, for competitive interactions with a continuous range of strategies we find that among-group selection may fail to produce socially optimal outcomes, but it can nonetheless produce second-best solutions that balance individual incentives to defect with the collective incentives for cooperation. We conclude by describing the broad applicability of multiscale evolutionary models to problems ranging from the production of diffusible metabolites in microbes to the management of common-pool resources in human societies.

Host heterogeneity and epistasis explain punctuated evolution of SARS-CoV-2.

Identifying drivers of viral diversity is key to understanding the evolutionary as well as epidemiological dynamics of the COVID-19 pandemic. Using rich viral genomic data sets, we show that periods of steadily rising diversity have been punctuated by sudden, enormous increases followed by similarly abrupt collapses of diversity. We introduce a mechanistic model of saltational evolution with epistasis and demonstrate that these features parsimoniously account for the observed temporal dynamics of inter-genomic diversity. Our results provide support for recent proposals that saltational evolution may be a signature feature of SARS-CoV-2, allowing the pathogen to more readily evolve highly transmissible variants. These findings lend theoretical support to a heightened awareness of biological contexts where increased diversification may occur. They also underline the power of pathogen genomics and other surveillance streams in clarifying the phylodynamics of emerging and endemic infections. In public health terms, our results further underline the importance of equitable distribution of up-to-date vaccines.

Social dilemmas of sociality due to beneficial and costly contagion.

Levels of sociality in nature vary widely. Some species are solitary; others live in family groups; some form complex multi-family societies. Increased levels of social interaction can allow for the spread of useful innovations and beneficial information, but can also facilitate the spread of harmful contagions, such as infectious diseases. It is natural to assume that these contagion processes shape the evolution of complex social systems, but an explicit account of the dynamics of sociality under selection pressure imposed by contagion remains elusive. We consider a model for the evolution of sociality strategies in the presence of both a beneficial and costly contagion. We study the dynamics of this model at three timescales: using a susceptible-infectious-susceptible (SIS) model to describe contagion spread for given sociality strategies, a replicator equation to study the changing fractions of two different levels of sociality, and an adaptive dynamics approach to study the long-time evolution of the population level of sociality. For a wide range of assumptions about the benefits and costs of infection, we identify a social dilemma: the evolutionarily-stable sociality strategy (ESS) is distinct from the collective optimum-the level of sociality that would be best for all individuals. In particular, the ESS level of social interaction is greater (respectively less) than the social optimum when the good contagion spreads more (respectively less) readily than the bad contagion. Our results shed light on how contagion shapes the evolution of social interaction, but reveals that evolution may not necessarily lead populations to social structures that are good for any or all.

Models with higher effective dimensions tend to produce more uncertain estimates.

Mathematical models are getting increasingly detailed to better predict phenomena or gain more accurate insights into the dynamics of a system of interest, even when there are no validation or training data available. Here, we show through ANOVA and statistical theory that this practice promotes fuzzier estimates because it generally increases the model's effective dimensions, i.e., the number of influential parameters and the weight of high-order interactions. By tracking the evolution of the effective dimensions and the output uncertainty at each model upgrade stage, modelers can better ponder whether the addition of detail truly matches the model's purpose and the quality of the data fed into it.

Coexistence in diverse communities with higher-order interactions.

A central assumption in most ecological models is that the interactions in a community operate only between pairs of species. However, two species may interactively affect the growth of a focal species. Although interactions among three or more species, called higher-order interactions, have the potential to modify our theoretical understanding of coexistence, ecologists lack clear expectations for how these interactions shape community structure. Here we analytically predict and numerically confirm how the variability and strength of higher-order interactions affect species coexistence. We found that as higher-order interaction strengths became more variable across species, fewer species could coexist, echoing the behavior of pairwise models. If interspecific higher-order interactions became too harmful relative to self-regulation, coexistence in diverse communities was destabilized, but coexistence was also lost when these interactions were too weak and mutualistic higher-order effects became prevalent. This behavior depended on the functional form of the interactions as the destabilizing effects of the mutualistic higher-order interactions were ameliorated when their strength saturated with species' densities. Last, we showed that more species-rich communities structured by higher-order interactions lose species more readily than their species-poor counterparts, generalizing classic results for community stability. Our work provides needed theoretical expectations for how higher-order interactions impact species coexistence in diverse communities.

Sociocultural determinants of global mask-wearing behavior.

Behavioral responses influence the trajectories of epidemics. During the COVID-19 pandemic, nonpharmaceutical interventions (NPIs) reduced pathogen transmission and mortality worldwide. However, despite the global pandemic threat, there was substantial cross-country variation in the adoption of protective behaviors that is not explained by disease prevalence alone. In particular, many countries show a pattern of slow initial mask adoption followed by sharp transitions to high acceptance rates. These patterns are characteristic of behaviors that depend on social norms or peer influence. We develop a game-theoretic model of mask wearing where the utility of wearing a mask depends on the perceived risk of infection, social norms, and mandates from formal institutions. In this model, increasing pathogen transmission or policy stringency can trigger social tipping points in collective mask wearing. We show that complex social dynamics can emerge from simple individual interactions and that sociocultural variables and local policies are important for recovering cross-country variation in the speed and breadth of mask adoption. These results have implications for public health policy and data collection.

Sharp thresholds limit the benefit of defector avoidance in cooperation on networks.

Consider a cooperation game on a spatial network of habitat patches, where players can relocate between patches if they judge the local conditions to be unfavorable. In time, the relocation events may lead to a homogeneous state where all patches harbor the same relative densities of cooperators and defectors, or they may lead to self-organized patterns, where some patches become safe havens that maintain an elevated cooperator density. Here we analyze the transition between these states mathematically. We show that safe havens form once a certain threshold in connectivity is crossed. This threshold can be analytically linked to the structure of the patch network and specifically to certain network motifs. Surprisingly, a forgiving defector avoidance strategy may be most favorable for cooperators. Our results demonstrate that the analysis of cooperation games in ecological metacommunity models is mathematically tractable and has the potential to link topics such as macroecological patterns, behavioral evolution, and network topology.

A PDE Model for Protocell Evolution and the Origin of Chromosomes via Multilevel Selection.

The evolution of complex cellular life involved two major transitions: the encapsulation of self-replicating genetic entities into cellular units and the aggregation of individual genes into a collectively replicating genome. In this paper, we formulate a minimal model of the evolution of proto-chromosomes within protocells. We model a simple protocell composed of two types of genes: a "fast gene" with an advantage for gene-level self-replication and a "slow gene" that replicates more slowly at the gene level, but which confers an advantage for protocell-level reproduction. Protocell-level replication capacity depends on cellular composition of fast and slow genes. We use a partial differential equation to describe how the composition of genes within protocells evolves over time under within-cell and between-cell competition, considering an infinite population of protocells that each contain infinitely many genes. We find that the gene-level advantage of fast replicators casts a long shadow on the multilevel dynamics of protocell evolution: no level of between-protocell competition can produce coexistence of the fast and slow replicators when the two genes are equally needed for protocell-level reproduction. By introducing a "dimer replicator" consisting of a linked pair of the slow and fast genes, we show analytically that coexistence between the two genes can be promoted in pairwise multilevel competition between fast and dimer replicators, and provide numerical evidence for coexistence in trimorphic competition between fast, slow, and dimer replicators. Our results suggest that dimerization, or the formation of a simple chromosome-like dimer replicator, can help to overcome the shadow of lower-level selection and work in concert with deterministic multilevel selection in protocells featuring high gene copy number to allow for the coexistence of two genes that are complementary at the protocell level but compete at the level of individual gene-level replication. These results for the PDE model complement existing results on the benefits of dimerization in the case of low genetic copy number, for which it has been shown that genetic linkage can help to overcome the stochastic loss of necessary genetic templates.

Immune Heterogeneity and Epistasis Explain Punctuated Evolution of SARS-CoV-2.

Identifying drivers of viral diversity is key to understanding the evolutionary as well as epidemiological dynamics of the COVID-19 pandemic. Using rich viral genomic data sets, we show that periods of steadily rising diversity have been punctuated by sudden, enormous increases followed by similarly abrupt collapses of diversity. We introduce a mechanistic model of saltational evolution with epistasis and demonstrate that these features parsimoniously account for the observed temporal dynamics of inter-genomic diversity. Our results provide support for recent proposals that saltational evolution may be a signature feature of SARS-CoV-2, allowing the pathogen to more readily evolve highly transmissible variants. These findings lend theoretical support to a heightened awareness of biological contexts where increased diversification may occur. They also underline the power of pathogen genomics and other surveillance streams in clarifying the phylodynamics of emerging and endemic infections. In public health terms, our results further underline the importance of equitable distribution of up-to-date vaccines.